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Abstract 

River streamflows are excellent climatic indicators since they integrate precipita- 
tion over large areas. Here we follow up on our previous study of the influence of 
solar activity on the flow of the Parana River, in South America. We find that the 
unusual minimum of solar activity in recent years have a correlation on very low 
levels in the Parana's flow, and we report historical evidence of low water levels 
during the Little Ice Age. We also study data for the streamflow of three other 
rivers (Colorado, San Juan and Atuel), and snow levels in the Andes. We obtained 
that, after eliminating the secular trends and smoothing out the solar cycle, there 
is a strong positive correlation between the residuals of both the Sunspot Number 
and the streamflows, as we obtained for the Parana. Both results put together imply 
that higher solar activity corresponds to larger precipitation, both in summer and 
in wintertime, not only in the large basin of the Parana, but also in the Andean 
region north of the limit with Patagonia. 
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1 1 Introduction 



2 Usually, studies focusing on the influence of solar activity on climate have con- 

3 centrated on Northern Hemisphere temperature or sea surface temperature. 

4 However, climate is a very complex system, involving many other important 

5 variables. Recently, several studies have focused in a different aspect of climate: 

6 atmospheric moisture and related quantities like, for example, precipitation. 
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Perhaps the most studied example is the Asian monsoon, where correlations 
between solar activi ty and precipita tion have been found in several time 
scales. For example, iNeff et al\ (120011 ) found strong coheren ce between solar 



variab ility and the monsoon in Oman between 9 and 6 kyr ago. lAgnihotri et al. 



( 120021 ) found that Indian monsoon intensity followed the sola r irradiance vari 



ability on centennial time scales during the last millennium. iFleitmann et al. 



( 120031 ) studied Holocene forcing of the Indian monsoon, and found that inter- 
vals of weak (strong) solar activity correlates with periods of low (hig h) mon- 
soon precipitation. On shorter time scales, iMehta and Laul ( 1997| ). found 
that, at decadal-multidecadal time scales, the correlation between the El Nino 
3 index and the monsoon rainfall is stronger when solar irradiance is above 
normal and viceversa. Correlations between solar activity and Indian mon- 
soon i n dec adal tini e scale s were also found by iBhattacharyya and Narasimha 
( I2OO5I ) and iKoderal (l2004j ). among others. 



Wang et al\ (|2005[ ) studied the monsoon in southern China over the past 9000 
years, and found that higher solar irradiance corresponds to stronger mon- 
soon. They proposed that the monsoon responds almost immediately to solar 
changes by rapid atmospheric responses to solar forcing. 
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All these studies reported a positive correlation, where periods of higher solar 



activit y correspond to periods of larger precipitation. In contrast. iHong et al. 



(i200ll ) studied a 6000-year record of drought and precipitation in northeast- 
ern China, and found that most of the dry periods agree well with stronger 
solar activity and viceversa. In the American continent, droughts in the Yu- 
catan Peninsula have been associated with periods of h igh solar activity an d 



have even been proposed to explain the Mayan decline (iHodell et alV l200l[ ). 



In th e same sense, st udies based 1 1 the w ater level of Lakes Naivasha ( iVerschuren et al 
2000) and Victoria ( IStager et ali l2005l ) in East Africa, report severe droughts 
during phases of high solar activity and increased precipitation during peri- 
ods of low solar irradiation. To explain these differences it has been proposed 
that increased solar irradiation causes more evaporation in equatorial regions, 
enhancing the n et transport of moistur e flux to the Indian sub-continent via 
monsoon winds (lAgnihotri et all |2002[ ) . 



However, these relationships seem to have reversed sign around 200 years ago, 
as severe droughts developed over much of tropical A f rica d uring the Dal- 



ton sunspot r ninirn um, ca. AD 1800-1820 ( IStager et ali l2005l ). Furthermore 



Stager et al\ ( 120071 ) studied recent water levels in Lake Victoria, and found 
that peaks in the ~ 11-year sunspot cycle were accompanied by water level 
maxima throughout the 20th century, due to the occurrence of positive rain- 
fall anomalies ~1 year before solar maxima. Similar patterns also occurred in 
at least five other East African lakes, indicating that these sunspot-rainfall 
relationships were broadly regional in scale. 
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A different approach was taken by lMauas and Flamenco! (120051 ) who proposed 
to study the streamffow of a large river, the Parana in southern South Amer- 
ica, as an indicator of precipitation. In fact, flows of continental-scale rivers 
are excellent climatic indicators since they integrate precipitation, infiltra- 
tions and evapotranspiration over large areas and smooth out local variations. 
Signals of solar activity h ave recently been found with spectral analysis tech- 
niques in the river Nile by iRuzmaikin et al\ ( 120061 ) , who found a low-frequency 
88- year variation presen t in solar variability and in the Nile records. Similarly, 
Zanchettin et al.\ ( 120081 ) found that the discharge of the Po river appear to 
be correlated with variations in solar activity, on decadal time scales. 



In Mauas et al. ( 20081 ) (hereinafter Paper I) we presented the results of our 
study of the Parana. We found that the streamfiow variability of the Parana 
river has three temporal components: on the secular scale, it is probably part 
of the global climatic change, which at least in this region of the world is 
related with more humid conditions; on the multidecadal time scale, we found 
a strong correlation with solar activity, as expressed by the Sunspot Number, 
and therefore probably with solar irradiance, with higher activity coincident 
with larger discharges; on the yearly time-scale, the dominant correlation is 
with El Nino. 



67 In the present paper we follow up on the study of the influence of solar activity 

68 on the flow of South American rivers. In Section 2 we expand in time the study 

69 of the multi-decadal component of the Parana's streamfiow, to include the 

70 most recent years, which have shown particularly low levels of solar activity. 

71 In Section 3 we study other South American rivers, to see whether the influence 

72 extends to other areas of the continent. Finally, in Section 4 we discuss the 

73 implications of our findings. 



74 2 The multidecadal component of the Parana's streamfiow 



75 The Parana is the fourth river of the world according to streamfiow (20 600 

76 m^/s), and the fifth according to drainage area (3 100 000 km^), which is the 

77 second largest in South America. Its origin is in the southernmost part of the 

78 Amazon forest, from where it flows south collecting water from territories in 

79 Brazil, Bolivia, Paraguay and Argentina. Its outlet is in the Plata River, a few 

80 kilometers north of the City of Buenos Aires. It flows through heavily popu- 

81 lated areas and it is navigated by overseas trade ships, unlike other rivers of 

82 similar size. For these reasons, its streamfiow has been measured continuously 

83 during the last century. 

84 As in Paper I, we analyze the streamfiow data measured daily since 1904, at 

85 a gauging station located in the city of Corrientes, 900 km north of the outlet 
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B6 of the Parana. Since the Parana's hydrological year goes from September to 

B7 August, with maximum streamflow in the Southern Hemisphere's summer 

88 months of January, February and March, our yearly series integrates the flow 

89 from September to August of the next year. 

90 In Paper I we found that in intermediate scales of decades, there is a strong 

91 correlation between the Parana's streamflow and solar activity, as expressed 

92 by the Sunspot Number (S'ArjZI, with larger solar activity corresponding to 

93 larger streamflow. We found a simila rly strong correlation with the irradiance 



95 We further explore this correlation in Fig. [H which is an update of Fig. 2 of 

96 Paper I, including 4 more years of data. To retain only the intermediate scale, 

97 we first computed the secular trends with a low-pass Fourier filter with a 50 

98 years cut-off, as shown in Fig. 1 of Paper I, which was substracted from the 

99 data. Then, we performed an 11-year-running mean to smooth out the solar 

100 cycle (for this reason, only data for the period 1909-2003 are shown). In this 

101 way, both high and low frequencies have been filtered out of the data in Fig. 

102 [H which only retain the variations in timescales between 11 and 50 years. 

103 When plotting together different quantities, the offset and the relative scales 

104 are free parameters which are usually arbitrarily introduced. To avoid these 

105 two artificial parameters, as a final step we have standardized the quantities 

106 by subtracting the mean and dividing by the standard deviation of each series 

107 shown, for the whole period 1909-2003. More details can be found in Paper I. 

108 The correlation between the Parana's streamflow and the Sunspot Number 

109 in Fig. [T]is quite remarkable. In fact, the correlation coefficient between both 
no series is R=0.78, significant to a 99% level. 

111 We point out that this correlation is found in the intermediate time scale. On 

112 longer timescales, both the Parana's discharge and solar activity are larger in 

113 the last decades than in the first ones of the 20th century, and these increases 

114 are not correlated (for a discussion, see Paper I). On the yearly timescale, 

115 the dominant factor influencing streamflow's variations is El Ninno (again, 

116 see Paper I for details). The results shown in Fig. [T] show that decades of 

117 larger discharge correspond to decades of higher activity, with these variations 

118 overimposed on the corresponding secular trends. 

119 It can be seen that the correlation is still found in the most recent years. 

120 In particular, in the period 1995-2003 both the mean Sunspot Number and 

121 the Parana's streamflow have decreased by similar amounts. In fact. Solar 

122 Cycle 23 was the weakest since the 1970s, and the onset of Solar Cycle 24 

123 was delayed by a minimum with the largest number of spotless days since the 



Available at ftp;//ftp.ngdc.noaa.gov/STP/SOLAR_DATA/SUNSPOT_NUMBERS. 
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Fig. 1. The detrended time series for the Parana's streamflow (full line) and the 
Sunspot Number (dashed line). The detrended series were obtained by subtracting 
from each data series the corresponding secular trend and were smoothed by an 
11-yr-running mean to eliminate the solar cycle. Both series were standardized by 
subtracting the mean and dividing by the standard deviation, to avoid introducing 
arbitrary free parameters. The Pearson's correlation coefficient is R=0.78. 

124 1910s. Morover, for the years 2008 and 2009 (2.9 and 3.1, respectively), 

125 have been the lowest since 1913. Similarly, the mean levels of the Parana were 

126 also the lowest since the 1970s. 

127 We have also tested the correlations between the Parana's discharge and the 

128 neutron count at Climax, Coloradc[l], which is a direct measure of the galactic 

129 cosmic rays (GCR) flux. Of course, since neutron count is correlated with 

130 sunspot numbers, we found a correlation with Parana's streamflow. However, 

131 the correlation with Sn is larger, pointing to a more direct correlation with 

^ Available at http;/ /www.env.sci.ibaraki.ac.jp/ftp/pub/WDCCR/STATIONS/climax. 
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The relationship between smaller solar activity and low Parana's discharge can 
also be found in historical records. For example, low discharges were reported 
during the period known as the Little Ice Age (LIA). In particular, a traveler 
of that period recalls in his diary that in the year 1752 the streamflow was 
so small that the river could not even be navigated by the ships of that time, 
which were less than 5 ft draft, to be compared with ships up to 18 ft draft 
that can navigate it at present as far north as Asuncion in Paraguay (jlriondd . 
Il999h . The fact that the LIA coincided with reduced pr ecipitation in this 



region has been found in different climatic records (e.g. iPiovano et al.\ 12009 
and references therein). It is well known that the LI A coincided with, and 
perhaps was caused by, low solar activity ( jEddyl . Il976l ). 



144 3 The Colorado river basin 



145 Here we study the streamflow of the Colorado river, and two of its tributaries, 

146 the San Juan and the Atuel rivers (see Fig. We also analyze snow levels, 

147 measured near the sources of the Colorado. 
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Fig. 2. Colorado hydrologic system. The rivers under study are marked in the figure: 
a. Colorado, b. San Juan, c. Atuel and d. the lower part of the Parana river. The 
stream flow (■) and snow (•) measuring stations are also indicated. 
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148 The Colorado river marks the north boundary of the Argentine Patagonia, 

149 separating it from the Pampas, to the northeast, and the Andean region of 

150 Cuyo, to the Northwest. Its origin is on the eastern slopes of the Andes Moun- 

151 tains, from where it flows southeast until it discharges in the Atlantic Ocean. 

152 The Atuel, which originates in the glacial Atuel Lake, at 3250 m above sea 

153 level in the Andes range, and the 500 km long San Juan river, join the Col- 

154 orado downstream of it's gauging station. Therefore, the data given by the 

155 three series are not directly related. 

156 In Table [T] we list the mean stream flow and drainage basin area of the Col- 

157 orado, Atuel and San Juan rivers. We also include the time interval of the 

158 stream flow records plotted in Fig. H] and the geographical coordinates of the 

159 gauging stations. 

160 Unlike the Parana, whose streamflow is directly related to precipitation, 

161 the regime of all these rivers is dominated by snow melting, and their stream- 

162 flows reflect precipitation accumulated during the winter, and melted during 

163 spring and summer. For this reason, the streamflows are largest during sum- 

164 mer, and the hydrological year for these rivers goes from July to June next 

165 year. This can be seen In Fig. |31 where we show the mean monthly flow of 

166 the Colorado. In the figure we separately plot the flow for years when the 

167 multidecadal component of the sunspot number shown in Fig. [T]is high (low), 

168 i.e. larger (smaller) than 0.5 a above (below) the mean value. It can be seen 

169 that during the decades with larger activity, the streamflow is larger from 

170 September to December, when most of the melting takes place, and remains 

171 almost constant during the rest of the year. 

172 To directly study the snow precipitation, we complete our data with mea- 

173 surements of the height of snow accumulated at Valle Hermoso (35° 15' S; 70° 

174 20' W), in the Andes at 2250 m above Sea level, close to the origin of the 

175 Colorado (see Fig. [2]), which were measured in situ at the end of the winter 

176 since 1952. In fact, the correlation between the streamflow of the Colorado and 

177 the snow height is very good, with a correlation coefficient R=0.87, significant 

178 to a 99% level. Correlation coefficients between the snow data and the Atuel 

179 and San Juan streamflows are R=0.76 and R=0.64, respectively. Since Valle 
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Fig. 3. Monthly Colorado's mean streamflow. The solid (dashed) curve was obtained 
considering only the years when the multidecadal component of Sn shown in Fig. 
[T]is high (low), i.e. larger (smaller) than 0.5 a above (below) the mean value. 

180 Hermoso is placed closer to the origin of the Colorado, and closer to the Atuel 

iBi than to the San Juan, this progressive reduction of the correlation is to be 

182 expected. 

183 In Fig. m we plot the yearly time series of the streamflow of the Colorado, San 

184 Juan and Atuel rivers, the snow height and the Sunspot number. We also show 

185 the variation in the longest scales, obtained with a low-pass Fourier filter, as 

186 we did for the Parana. However, since the length of the time series is not the 

187 same for every set of observations, we could not apply a uniform filter for all of 

188 them. In all cases, the cut-off was taken as half the length of the observations 

189 (33 years for the Colorado, 40 years for the Atuel, 50 years for the San Juan 

190 and the Sunspot Number and 28 years for the snow) 

191 In Fig. |5]we compare the multidecadal component of the streamflows with the 

192 corresponding series for the sunspot number. In each case, we smoothed out 

193 the solar cycle with an 11-year running mean, and we detrended the series by 

194 subtracting the long term component shown in Fig. HI Finally, we standard- 

195 ized the data by subtracting the mean and dividing by the standard deviation 

196 of each series shown for the period 1971-2000, suggested by the World Mete- 

197 orological Organization as standard reference. In the panel corresponding to 

198 the Colorado, we also include the snow height. 

199 It can be seen that in all cases the agreement is remarkable. The correlation 

200 coefficients are 0.59, 0.47, 0.67 and 0.69 for the Colorado, the snow level. 
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Fig. 4. Time series of the Colorado, San Juan and Atuel's streamflow, the snow 
height measured at Valle Hermoso and the Sunspot Number. The long-term trend 
of each series is marked with a heavy line. 

201 the San Juan and the Atuel, respectively, all significant to the 96-97% level. 

202 Therefore, also in these cases we found a relation between solar activity, on 

203 one hand, and the streamflow of these rivers and snow level, on the other, as 

204 we found for the Parana. 



Probably, the most important correlation is the one with snow level, 
and the correlations with the rivers' streamflow are indirect conse- 
quences of the variations in precipitations. We should point out that 
climate in this area is correlated with the conditions over the equatorial Pa- 
cific, as measured by El Nino. This correla tion was studied for t he Di amante 
River, also a tributary of the Colorado, by iBerri and Flamenco! (1l999[ l. 



211 In particular, the peaks in the snow level and the streamflows (see Fig. S]) in 
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212 the year 1982/3 coincide with a very strong El Niiio event, which caused a huge 

213 flood in the Parana's basin, as we discussed in Paper I. Correlation coeflicients 

214 between our data and el Nino 1+2 index in November (at the beginning of the 

215 austral Summer), are R=0.51, R=0.60, R=0.60, and R=0.61, for snow level 

216 and the Colorado, Atuel and San Juan rivers, respectively. 

217 Although all these rivers have maximum streamflow during Summer, there is a 

218 big difference, however, between the regimes of the Parana and the remaining 

219 rivers: for these ones, the important factor is the intensity of the precipitation 

220 occurring in the winter months, from June to August. For the Parana, what is 

221 most important is the level of the precipitation during the summer months. 

222 It is also worth noticing the sense of the relationship: here again, stronger 

223 activity coincides with larger precipitation. 



224 4 Discussion 

225 In this paper we analyzed the influence of solar activity in the streamflow 

226 of South American rivers of different regimes. First, we extended in time the 

227 study of the correlation between Sunspot Number and the Parana's streamflow 

228 we reported in Paper I. On one hand, we found that the unusual minimum 

229 of solar activity in recent years have a correlation on very low levels in the 

230 Parana's flow. On the other, we reported historical evidence of low water 

231 levels during the Little Ice Age. We also found that the correlation is stronger 

232 with sunspot number than with neutron count, which confirms that what is 

233 affecting climate is most probably solar irradiance, and not GCRs. 

234 The fact that the river's behaviour follows Sn through one more minimum 

235 strongly enhances the significance of the correlation and its predictive value. 

236 In particular, the low levels of activity expected for Solar Cycle 24 anticipate 

237 that the dry period in the Parana will continue well into the next decade. 

238 To study whether the solar infiuence extends to other areas of the continent, 

239 we analyzed the streamfiow of three South American rivers: the Colorado 

240 and two of its tributaries, the San Juan and Atuel rivers. We also used snow 

241 level from a station at the origin of the Colorado. We obtained that, after 

242 eliminating the secular trends and smoothing out the solar cycle, there is a 

243 strong correlation between the residuals of both the Sunspot Number and the 

244 streamfiows. In all cases, the correlation we found on multi-decadal time scales 

245 is positive, i.e., higher solar activity corresponds to larger snow accumulation 

246 and, therefore, to larger discharges of all these rivers, as we obtained for the 

247 Parana river. 

248 Therefore, both results put together imply that higher solar activity corre- 
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Fig. 5. The detrended streamflows (full 
lines) compared with Sunspot Number 
(dashed lines). In panel (a) the snow level 
is also shown (dotted line) . The data were 
smoothed with an 11-year running mean, 
detrended by substracting the long term 
component. All series were standardized 
by subtracting the mean and dividing 
by the standard deviation of each series 
shown for the period 1971-2000. 



249 spends to larger precipitation, not only in the large basin of the Parana, but 

250 also in the Andean region north of the limit with Patagonia. Furthermore, 

251 since streamfiow variability of rivers on central Chile are controlled by the 

252 same mechanisms that regulate the rivers studied in this paper, one might 

253 expect the same correlation to be found west of the Andes. 



Solar activity can affect precipitation through the position of the Inter Trop- 
ical Convergence Zone (I TCZ), which has been shown to corre late with varia- 
tions in solar insolation ( jPoore et ali l2004l : iHaug et all l200ll ). In fact, it has 
been proposed that a displacement southwards of the ITCZ would inc rease 
precipitation in southern tropical South America (INewton et a/.l . l2006l ) . We 
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259 point out that increased precipitation occur both in the Southern Hemi- 

260 sphere's summer when the ITCZ is over the equator, close to the origin of the 

261 Parana, and in wintertime, when the ITCZ displaces north, and precipita- 

262 tion increase further South. 
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